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A Semi-Empirical Model of the Infrared Emission from Galaxies 
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ABSTRACT 

We present a semi-empirical model for the infrared emission of dust around star-forming 
sites in galaxies. Our approach combines a simple model of radiative transfer in dust clouds 
with a state-of-the-art model of the microscopic optical properties of dust grains pioneered 
by Draine & Li. In combination with the STARBURST99 stellar spectral synthesis package, 
this framework is able to produce synthetic spectra for galaxies which extend from the Lyman 
limit through to the far-infrared. We use it to probe how model galaxy spectra depend upon 
the physical characteristics of their dust grain populations, and upon the energy sources which 
heat that dust. We compare the predictions of our model with the 8- and 24-|xm luminosities 
of sources in the Spitzer First Look Survey, and conclude by using the models to analyse the 
relative merits of various colour diagnostics in distinguishing systems out to a redshift of 2 
with ongoing star formation from those with only old stellar populations. 



Key words: infrared: galaxies - infrared: ISM 
ies: starburst - galaxies: high-redshift 



dust, extinction - radiative transfer - galax- 



1 INTRODUCTION 

It is increasingly clear that the infrared spectra of galaxies hold vi- 
tal clues concerning galaxy energetics. Observations of our nearest 
neighbours tell us that around 60 per ce nt of their star format ion is 
obscured by dust at visible wavelengths ( Fakeuchi et al. 2006| ), and 
the cosmic infrared background (GIB) indicates by its brightness 
that this is_alsoJrue_o£_sources_a^_co^ 

& D wek [200 Furthermore, deep sub-millimetre surveys (Smal l 



et al. |1997| ; [Hughes et alj|1998| ; |Smail et alj|1998| ; ^lain etal]|l999| ) 
have revealed a larg e population of ul tra-luminous infrared galaxies 
(ULIRGs) at z ~ 2 ( Blain et al. 2002 1 - clearly very dusty systems. 



and, if not harbour ing active nuclei, very actively star-forming also 
{^mail et alj[l997[). Given the degree of optical extinction in these 



dominates, m arked by a series of broad emission features (see, e.g., 
Draine|2003 , and references therein), the catalogue of which has re- 
centl y been greatly expanded by Spitzer ([Praine & Li[|2007[; S mith 
et al. [200"/[ ). These are attributed to polycyclic aromatic hydrocar- 
bon (PAH) molecules, and so we shall refer to them as 'PAH fea- 
tures'Q Finally, in the rest-frame near-infrared (NIR) - here taken 
to refer to A < 4 ii.m - stellar emission dominates. 

This wealth of available information has motivated many stud- 
ies which have sought to provide a framework in which this emis- 
sion can be interpreted. Some of these take an empirical approach, 
matching unresolved sou rces to template spectra derived fro m a 
vari ety of local galaxies ([Rowan-Robinson & Crawforj[l989|; Xu 



systems ( 5winbank et al. 2004 ), it is apparent that UV/optical stud- 



ies of the cosmic star formation history are subject to substantial 
incompleteness. 

The advent of the Spitzer Space Telescope has allowed a much 
greater understanding of the sources which comprise the GIB. Its 
resolution and sensitivity has allowed more than two million in- 
frared galaxies to be resolved in 49deg^ of sk y by the 5piYzer Wide 



et al. [2004 [Rowan-Robinson et al.|2004 [2005| ). These yield fast di- 
agnostics which are readily applicable to large numbers of sources. 
But they provide little information about the physical processes 
which fundamentally shape the spectra. 

Others have sought to develop semi-empirical models of in- 
frared spectra, co nsidering the propagation of radiation th rough 
dusty media (e.g. [Silva et al. || l994 l E fstathiou et al.j [20001: Tak- 



Area Infrared Extragalactic (SWIRE) Survey ( Lonsdale et al 



2003 



2004). Moreover, Spitzer s wavelength coverage, 3.6-160 nm, en- 
compasses three emission regimes in the spectra of normal galax- 
ies, each yielding information complementary to the others. In the 
rest-frame far-infrared (FIR) - here taken to extend from around 
30 to 300 \im - thermal emission from large dust grains dominates. 
In the rest-frame mid-infrared (MIR) - here taken to extend from 
around 4 to 30 \im ~ emission from transiently-heated dust grains 

* Email: dcf21@mrao.cam.ac.uk 



agi et al. 2003 ; Efstathiou & Rowan-Robinson 1995 ). This task can 
be split into two parts: modelling the optical properties of individ- 
ual dust grains, and modelling the large-scale transport of radia- 
tion throughsome^ealistic^ 
et al . ^006f '' ^iovai7'er'al] ' p006a| ]b^^ 

consideration of the radiative transport, including a treatment not 
only of absorption, but also of elastic photon scatterin g, which be- 
comes the principal source of complexity. Others (e.g. \A & Draine 



' In the literature, they are also commonly referred to as 'Aromatic Infrared 
Bands' (AIBs), or, historically, as 'Unidentified Infrared Bands' (UIBs) 
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200l| ; |Draine&Li|2007| ) tiave developed highly sophisticated mod- 
els of the optical properties of individual dust grains, but simplify 
radiative transfer. Essentially a computational trade-off exists: stud- 
ies which include a detailed treatment of radiative transfer involve 
many evaluations of dust emissivities in differing environments. 
This is prohibitively time consuming with state-of-the-art models 
of the microscopic optical properties of individual dust grains and 
so simpler alternatives must be sought. 

In this paper, we use a state-o f -the-a rt mo del for the ind i vidua l 
dust grains - that of Draine & Li (2001) and Li & Draine (2001) 



Ray Y 



- and then seek to make a minimal set of simplifications in our 
treatment of the radiative transfer problem such that it becomes 
computationally viable. 

In Section □ we outline the dust geometries which we consider. 
In Section |3| we go on to present our modelling of the radiative 
transfer, and in Section ^ we construct UV-visible radiation fields 
appropriate for dust-heating by star-forming galaxies. In Sections ^ 
and eI we draw upon models of dust grain populations and their 
microscopic optical properties from the literature. In Section ^ we 
conclude the development of our model with a simple framework 
for modelling the evolution of the metallicity and mass of gas in 
passively-evolving star-forming gas clouds. In Sections | 8| a nd ^ we 
present the basic predictions of our model. In Section lid we use 
the model to predict 8- and 24-jj.m luminosities for a sample of 
galaxies and compare to observations. Finally, in Section MJ, we 
develop a simple model for the evolving colours of early- and late- 
type galaxies. 

Wherever required, we assume a flat Acdm cosmology with 
Hq = 72kms"' Mpc"' and JIa = 0.7. 



2 THE MODEL GEOMETRY 

In this paper we consider two basic geometries for the spatial distri- 
bution of dust and the source of illumination. The first geometry is 
a shell of dust grains - which we shall term a 'circumnuclear' grain 
population ~ surrounding a point-like heating source, representing 
dust around a star-forming region. The second is a uniform distri- 
bution of dust grains - which we shall term a 'diffuse' grain pop- 
ulation - spread throughout a diffuse inter-stellar medium (ISM) 
within which the radiation field is assumed spatially uniform. In 
a future paper, we shall use them as components of a composite 
model for star-forming galaxies. 

The geometry adopted for our diffuse populations is the sim- 
pler. We assume the dust-bearing ISM to be optically thin at in- 
frared wavelengths, such that the re-absorption of dust emission 
can be neglected. We further assume the UV-visible heating radia- 
tion field within it to be spatially uniform, hence all grains of any 
given size and composition have the same emissivity. Under these 
conditions, total dust emission is directly proportional to the num- 
ber of grains present, and we therefore scale all quantities per unit 
volume of ISM. 

The circumnuclear geometry is illustrated in Figure [l]. A 
point-like heating source lies at the centre of a spherical shell of 
dust, of inner radius and outer radius rj . Within the shell, we 
trace the density of the medium via the number density of hydrogen 
nuclei, assumed to be spherically symmetric and denoted nii{r), 
such that the column density A^c of hydrogen nuclei along a line of 
sight passing through a dust shell to its nucleus is given by: 

Ac= rnn{r)dr. (I) 

Astrophysically, the cavity at r < ro might correspond to an H II 



Dust Shell 

H number density nn{7') 
Grain number density ngr(r) 




Figure 1. Dust in a spherical shell of dust around a single isotropic heating 
source. 



region, essentially devoid of dust due to the sublimation of grains 
by energetic photons. 



3 RADIATIVE TRANSFER 

We adopt a highly-simplified treatment of the radiative transfer 
which we argue is sufficiently accurate for a wide range of prob- 
lems. 



3.1 The circumnuclear geometry 

The evolution of the surface brightness ly along ray Y in Figure |l| 
is governed by the time-independent equation of radiative transfer: 



dlv 
d7 



^ -Cv,ext'!H('')/v+ev('')nH('') + 



mir) 
4n 



Cv.sca(0')dii', 

(2) 

where the first term on the right-hand side describes the absorption 
of radiation by dust, the second dust emission, and the third the 
scattering of photons into the ray. The integral in the third term is 
over solid angle Q.' alX; 0' is the scattering angle between dfi' and 
the direction of Y. 

As mentioned above, we use «h('') to parameterise the varying 
spatial density of material. Cy.abs and Cy.sca are the cross sections 
to absorption and scattering respectively, expressed per hydrogen 
atom, and averaged as described below over the compositions and 
sizes of particles within the grain population. £v{r) is the emissiv- 
ity of the grain population - the power emitted per unit frequency 
into unit solid angle, normalised in the same way as for the cross 
sections, s measures distance along the ray. 

The extinction cross section, Cy.ext^ is the sum of the absorp- 
tion and scattering cross sections: 

Cv.ext = Cv.abs +Cv,sca- (3) 

The averaging of the quantities above over grains of varying com- 
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positions and sizes is performed as follows: 



1 d4,(fl) 



«j] da 



da, 



da. 



(4) 
(5) 



where j e {ext, abs, sea}, i denotes a population of grains of given 
composition, and «gr(a) denotes the spatial number density of 
grains of each composition with radius smaller than a. 

The first simplification that we make is to neglect scattering, 
such that Cv,sca = 0. We discuss the validity of this assumption 
in the following section. Integrating the remaining terms of Equa- 
tion (^) over all rays passing through the point X iri Figure |l[ we 



obtain ( Z^handresekhai 



1960 



Rowan-Robinson 



1980| ): 



1 

'{'') = - lydfi, fl = cosO. 

2 7-1 

1 /•' 



where: 



and: 



(6) 



(7) 



(8) 



with denoting the angle made between each ray and the radial 
direction, as shown in Figure [l| 

Geometrically, Jy (r) may be visualised as the average surface 
brightness along all rays passing through X, averaged over 47t stera- 
dians. Hy[r) may similarly be visualised as the average projected 
onto the radial directio n. 

Rowan-Robinson (1980) introduced what has become a 



widely-used decomposition of the surface brightness ly along each 
ray into three components, depending upon where photons last in- 
teracted with matter: 



'V — ly 



(9) 



where /y'' is radiation from the central heat source which has not 
been absorbed by dust, /y^' is radiation emitted by dust and I^^ is 
scattered radiation - which we have already neglected. While we 
do not use this in our mathematical treatment of Equation (^, it is 
useful in our discussion presently. 

The relationship between Jy (r) and Hy [r] encodes the angular 
distribution of the radiation flux passing through X. In the limiting 
case of a radiation field propagating exclusively in the radial direc- 
tion, Hv{r) = Jv{r). In the opposite limit of an isotropic radiation 
field, Hy{r) = i/2jy{r). Given a point heat source, as in Figure [l[ 

the former limit is applicable to the component /y''; the radiation 

(2) 

field emanating from the heat source is purely radial. For ly , how- 
ever, Hy{r) < Jy{r). 

Our second, and final, assumption, is that Jv{r) = Hy{r) in 
Equation (^. For UV-visible wavelengths, this assumption holds 
because the heating radiation field is expected to dominate over 
dust emission at these wavelengths, and so ly ~ /y' ' . In the infrared, 
where ly is expected to dominate ly, Jy{r) is under-predicted, but 
the assumption continues to hold if £y{r) 2> Cy_gi,sJv(r), that is 
to say, if the dust emission from the shell is not appreciably re- 
absorbed. Geometrically, this assumption is equivalent to assuming 
that dust emission is beamed along the outward radial direction. 

Finally, we note that the net outward flux Fy(r) of radiation 



through the sphere of constant radius passing through X is related 
to Hyij) via: 

Fy{r)=A7tHy{r). (10) 
Equation can thus be re-written: 

1^ ('-Vy(r)) = nn{r) {AnEy{r) -Cy,,k,Fy{r)) . (II) 

This equation is integrated numerically from the inner to the outer 
radius. 



3.2 Assumptions made in the circumnuclear geometry 

The assumptions made in the previous section - i.e. the neglect of 
scattering and the radial beaming of dust emission - will have neg- 
ligible effect upon the predictions of our model for column den- 
sities of dust which are optically thin at all wavelengths, i.e. for 
Nc < lO^-' Hm^^. For column densities of dust which are optically 
thick in the UV, but not in the infrared, our neglect of scattering will 
lead us to under-estimate the path lengths of UV/optical photons 
through the dust shell by a factor of 1-2, and so to under-estimate 
the absorption of UV/optical radiation by a similar factor. Since this 
effect is essentially the same as that of reducing the column density 
of dust, the effect when using these models to fit the spectral shape 
of real sources will be that we will over-estimate the dust masses 
of these objects. 

Our assumption that dust emission is beamed radially out- 
wards only begins to fail for dust shells with higher column den- 
sities still, when they become optically thick even at infrared wave- 
lengths. As the dust emission is assumed to take the shortest path 
out of the dust shell, we will under-estimate its re-absorption in 
these optically thick cases. In practice, this effect becomes signif- 
icant for dust column densities > 10^''Hm^^, as will be shown in 
Figure pi] 



Kowan-Kobmson (1980) studied the effects of a similar set of 



assumptions in their calculation of radiative transport in hot-centred 
star-forming clouds, and for the range of systems they consider, 
they report errors of around 10 per cent. 

In addition to the two assumptions just discussed, it is also 
apparent that the adopted geometry is simplistic, having only one 
single heat source. We note, however, that this geometry is observa- 
tionally indistinguishable from an ensemble of A' smaller circum- 
nuclear geometries, each heated by its own central heat source with 
luminosity scaled by a factor I /N with respect to the single shell, 
and each containing a dust shell with inner and outer radii scaled 
by factor I /\/N and dust density distribution n'{r) scaled according 
to: 



n(r) = VNn{r\/N) 



(12) 



where n{r) is the density distribution of the single shell. 

The column density of dust around each heat source in this lat- 
ter ensemble is the same as that in the former single circumnuclear 
shell; the radiation field incident upon grains on the inner edge of 
each dust shell is the same; the total mass of dust in the two cases is 
the same; and to a remote observer, the total solid angle subtended 
by the dust in the two cases is the same. In summary, although 
our circumnuclear model is nominally of dust around a single heat 
source, it is also a good model of systems where that luminosity 
production is distributed between several discrete sources. 
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3.3 The diffuse geometry 

Our treatment of radiative transfer in diffuse dust grain populations 
is simpler than the above. Instead of having a heat source of lu- 
minosity Lv> we have an interstellar radiation field (ISRF), whose 
energy density iiy'^ we normalise to XXo^ where Xo is that of the 
solar neighbourhood interstellar radiation field less the cosmic mi- 
crowave background (CMB)^ and is a free parameter. We adopt 
Xo = 7.46 X 10^'"* Jm^- ^, derived from integ ration of the ISRF of 



Mathis et al. 



(1983) and Mezger et al, (1982). 



The total luminosity emerging from the model can then be 
written: 



: 47t£v riiiV - 



(13) 



where V is the volume of the dust-bearing ISM, and A its surface 
area through which the interstellar radiation field leaves the galaxy. 



4 THE HEATING RADIATION FIELD 

In this paper, we consider models for the heating radiation field 
due to star formation. We generate these using Ve rsion 5.1 of the 



STARBURST99 Stellar spectral syn thesis package (Leitherer et al 



1999; Vazquez & Leitherei 2005). This package offers the facil 



ity to model stellar populations with two classes of star formation 
histories (SFHs). The first is ongoing star formation, proceeding at 
constant star-formation rate (SFR) "P, which began at some time 
t previously; we hereafter term these models 'continuous' SFHs. 
The second is a delta-function SFH, representing an instantaneous 
burst of star formation at some time t previously, of total mass m; 
we term this the 'instantaneous' SFH. The stellar populations mod- 
elled for the instantaneous SFHs may be referred to as single stellar 
populations (SSPs), which is to say that all of the stars within them 
are coeval. 

STARBURST99 models both stellar emission and also nebu- 
lar continuum. For the former it uses the stellar evolution tracks of 



the Padova group (Fagotto et al. 1994), with the addition of tracks 



for thermally pulsing asymptotic giant branch (TP-AGB) stars to 
improve the accuracy of the modelling of low and intermediate 
mass stars .r For the nebular continuum, the emission coefficients 



n^( 



of iFerlandl (198C) are used. This is the s ource of a problem in 
STARBURST99Q the data of ^erlanc ( 198C ) do not extend beyond 
4.5 (xm, however they are extrapolated to 160ii.m, yielding large, 
unphysical, infrared luminosities. For X ^ 4^m, we use a more 
physical extrapolation, taking Ly = v^" '. 

In addition to the continuous and instantaneous SFHs mod- 
elled by STARBURST99, we have also considered arbitrary SFHs, 
modelled by convolving the luminosity, Lggp (r) , of a 10^ M0 single 



^ We neglect the CMB in this nomialisation because, in contrast to the 
starlight component of the ISRF, it would make no sense to enhance it by 
a factor X- It should be noted that the CMB is also absent from all models 
presented in this paper. 

^ It should be noted that this is a significant departure from previous ver- 
sions of STARBURST99, which used the stellar evolution tracks of the 
Geneva group, and did not model low mass stars, introducing serious er- 
rors in the modelling of old stellar populations 



stellar population of age / with our star formation history: 



10* Mr: 



Lssp{t-t')dt'. 



(14) 



5 THE DUST MODEL 

Whilst there exist fairly tight observational constraints on the com- 
position and size distribution of dust grains in the Milky Way, rel - 
atively little is known about those in other galaxies ( Draine 2003 ). 
In this paper, we therefore base our dust grain population upon that 
inferred from o bservation of o u r own galaxy. 

Following Li & Draine (2001, hereafter, LDOl), we con- 
sider binary populations of 'carbonaceous' and silicate grains. The 
former sub-population includes both PAH molecules and larger 
graphitic grains; the optical properties of these grains exhibit a 
smooth transition with grain radius, centred around a radius of . 
The absorption cross sections of carbonaceous grains of radius a is 
taken to be: 

C-^bs(«) = ^PAH{a)Cl^SM + [1 - ^pah(«)]C^;L(«), (15) 
where the weighting parameter i^pah('^) is given by: 

^PAtiia) = (1-%-a) X min [l,(a^/a)^] , (16) 

and ^gia = 0.01 sets even the smallest PAH molecules to exhibit 1 
per cent of the continuum absorption of graphitic grains. The tran- 
sition radius, , is set by default to 50 A. 

For the size-distribution of grains in ea ch of these populations, 
we fo llow the parametric forms used by Weingartner & Draine 
(2001); for the carbonaceous grains, we use: 



1, 

exp 



dn|f 
da 



■-D{a)- 



a \at,e 



>(fl,j8g,a,,g) (17) 



3.5A < a < 



|- [(a-at,g)/ac,g]^|, flt.f 



«t,g 



< a 



where D{a) represents two log-normal peaks: 



D(a) 



-exp 



ln(a/ao,,) 



(18) 



which were introduced by LDOI to reproduce the mid-infrared lu- 
minosities observed by ISO, the Infrared Telescope in Space (IRTS) 
and by IRAS at 60 (xm. The term F{a,l5g, at g ) provides curvature: 



F{a,Pg,at) 



I + Pa/ at 
(l-pa/ai)- 



j8>0 
j8<0, 



(19) 



Z is the mass ratio of metals, Zr: 



0.02 is the solar mass ratio 



of metalSj_andal^^ther_s^mbols^^ & 
DraineT ^OoT^ I'Thrnormahsatioirc^^ by: 



Bi 



3 exp (-4.50^ 



pal.O 



l+erf|l|- 



ln(^£io.;/3.5A) 



(20) 



See notice by Hunt, October 30, 2006, in the Knowledge Base of the where niQ = 1.99 X 10 ^"kg is the mass of a carbon atom, p 



STARBURST99 website. 



2.24 X lO^kgm ^ is the density of graphite, ao.l = 3.5 A and 
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Parameter 


Value 




-1.54 


ft 


-0.165 


at,g 


0.0107 urn 


«c.E 


0.428 |jm 


Cg 


9.99 X 10- 


0!s 


-2.21 


ft 


0.300 


Ot,s 


0.164 urn 


Q 


1.00 X 10- 


be A 


0.75foc 


he. 2 


0.25foc 


he 


6.0 X 10-5 



Ta ble 1. The parameters of the default size distribution which we adopt - 
the Weingartner & Draine ( 2001 ) preferred distribution for i?v =3.1 Milky 



Way sight lines. 



ao.2 = 30 A are the wavelengths of the centres of the two log-normal 
peaks, and a = 0.4. 

For the silicate grains, we use: 



1 

1, 

exp 



da 



at.s 



F(a;j3s,at,s 



(21) 



|- [(fl-at.s)/ac,s]"'| : 



3.5A < a < at,s 

flt.s < CI 



By default, we s et the parameters of th e se siz e distributions 
to those preferred by Weingartner & Draine (ZOOl) for Ry = 3.1 
Milky Way sight lines, as given in Table |l|. The effect of using in- 
stead the size distributions preferred by those authors for Ry = 4.0 
and R\i =5.5 Milky Way sight lines will be discussed in Section]^ 

The appearance o f Z in Equation (|lj) is a departure from 
Weingartner & Draine (2001 ), who only considered Galactic envi- 
ronments. Thus, their size distributions are normalised to a dust-to- 
gas-mass ratio appropriate for solar metallicity environments. As 
the variation in this ratio with Z is quite poorly understood, we 
make the assumption that it is linearly proportional to Z, which is 
implicit in our renormalisation above. 

These distributions are plotted in Figure |^ To calculate the 
absorption cross sections of silicate and graphitic grains, we fol- 
low LDOl^nd^s^^diekcta^^nctims^ & 



Lee |1984| ) 'anT^ie'Theor^7seere!g!7 pohre^^ |l998| ) to 

estimate the absorption cross sections of spherical particles of ra- 
dius a. The treatment of graphitic grains is slightly complicated by 
the anisotropy of graphite's dielectric function. We follow LDOl in 
calculating an a veraged absorption cross section using the 'i/3-% 
approximation' ( Draine & Malhotn] 1993). 

For the PAH molecules, LDOl give algebraic fits to terrestrial 
labor atory measurements of for neutral and ionised sam- 

ples. Draine & Li ( 2007| , hereafter, DL07) revise these in the light 
of new near-infrared data ( Mattioda et al, 2005b I, and in order to fit 



the high-fidelity spectra of nearby star-forming galaxies observed 
by the S pitzer Infrared Nearby Galaxies Survey (SINGS; K enni- 
cutt et al.'^OOSp^rqjecr'^eTmpIe^^ by 
both LDOl and DL07, which are shown in Figures |(a) and |(b) for 
neutral and ionised grains respectively, both of radius 5 A. In the re- 
mainder of this paper, we use the DL07 cross sections throughout, 
except in Section llQ. In both cases, we average the cross sections 



100 



10 



i« 100 



10 



Carbonaceous grains 




10" 



10-'^ 
a I m 



10" 



10" 



Figure 2. The adopted size distributions for silicate (top) and carbonaceous 
(bottom) grains, assuming a solar metallicity environment. Unit areas under 
each distribution represent unit masses of grain material. 



of neutral and ionised PAH molecules with a weighting parameter 
/ describing the ionisation fraction. We t ake this to have a de fault 
value of 80 per cent, matching that which 



Draine & Li 



(2001) find 



in their model fits to Galactic photo-dissociation regions. The re- 
sulting absorption cross section C^J^l^^^ [a) is shown for a range of 
grain radii in Figure ^(c), and the absorption cross section C^'l,(,j. (a) 
of the silicate grain population Figure ^d). 



6 MODELLING THE EMISSIVITY OF THE DUST 

In this section, we outline how we model the emissivities ej, [a, r) 
of dust grains as a function of the energy density Ey [r) of radiation 
to which they are subjected. 



6.1 Transiently-heated grains 

Modelling emission from transiently-heated grains requires the cal- 
culation of the time-averaged probability distributions P{E) for 
their internal energies. Exact treatment of this problem would re- 
quire knowledge of all of their vibrational energy le vels and tran 



sition probabilities. Our simplified analysis follows Draine & Li 
(2001). We use Debye models for the normal modes of the C/Si 
skeletons o f PAH and sil i cate p articles, with Debye temperatures 
as used by Draine & Li (2001). We use Einstein models for the 
stretching, in-plane bending, and out-of-plane bending modes of 
the peripheral C-H bonds of PAH molecules: we assume the modes 
of each bond to be quantum harm onic oscillators with the same fun- 
damental frequencies, as given in Draine & Li (2001 ). 



To reduce the resulting mode spectra to a computationally 
tractab le number of energy states, we follow G uhathakurta & 
Draine ( 1989| ) and praine & Li ( ^001 ) in dividing them into A'bin 
bins (for our choice of bins, see Appendix A), with mean energies 



6 




Figure 3. The adopted absorption cross sections. Panels (a) and (b) show those for neutral and ionised PAH grains respectively, as given by LDOl and DL()7 
for grains of radius 5 A. A discussio n of the near-infrared f eature introduced into the cross sections of ionised PAHs by DL07 at 1.05 |jm, and the negative 
feature at 1.905 |jm, can be found in Mattioda et al. (2005a). Panels (c) and (d) show those for carbonaceous and silicate grains respectively for a variety of 
grain radii a, assuming a PAH ionisation fraction / = 0.8. Each trace is normalised with respect to the classical grain cross section of na- . 



Ui, widths AUi, and time-averaged occupation probabilities Pi. We 
denote as Tj,- the transition rate between bins i and j. The time 
evolution of Pj is then given by: 

dP. 



constramt: 



(22) 



The time-averaged steady-state probability distribution which we 
seek is that to which the above converges over time, and for which 
dPi/dt = 0. 

The elements of Ty; with j > i describe the upward transitions 
of grains that result fr om photon ab s orptio n; we model these using 
equations (15-25) of Draine & Li ( 2001 ); in Appendix B we re- 
produce these relations and describe a numerical optimisation that 
we use in their calculation. The elements with j < i describe the 
radiative cooling of grains; he re we use the 'the rmal continuous' 
approximation (equation 41 of Draine & Li 2001, reproduced here 
as Equation |b4| ), which models the cooling of grains as a continu- 
ous process, where each state i only makes downward transitions to 
the adjacent state i — 1. This allows much faster solution of Equa- 
tion E3 to find Pj. 

The diagonal terms are chosen ( [Praine & Li 2001 ) so that: 



hence Equation (E2|) can be re-written in the form: 



APi ^" 

7=0 



At 



(23) 



(24) 



These equations are solved, subject to the additional normalisation 



(25) 



using the method of Guhathakurta & Draine (1989). Given the vec- 



tor Pi, we calculate the time-averaged emissivity of each grain us- 
ing the thermal approximati on, under wh i ch ey {a,r) can be calcu- 
lated using equation (56) of Draine & Li ( 2001 ): 



2hv^ 



' sxp{hv/k6i) — 1 



(26) 



where 0; is t h e char acteristic temperature of bin as defined in 
Draine & Li (2001), the sum is over all bins / whose central 



energies are greater than hv, and we have neglected the factor 
(1 + X^ue /SJt) shown by those authors; this represents stimulated 
emission and may straightforwardly be shown to be negligible in 
all of the models presented in this paper. 



6.2 Large grains 

For sufficiently large grains, the approach outlined above becomes 
inefficient. Their internal energies become much larger than the en- 
ergies of the photons they absorb, and so their temperature fluctua- 
tions are not significant. Their internal e nergy probability d istribu- 
tions P{E) tend towards delta functions ( Li & Draine 2001 ). In this 
limit, we can model the energetics and emission of these grains by 
numerically solving the equation of radiative balance to find their 
equilibrium temperatures: 



/ Cv.abs(a)c"vdv= / 4;rCv,abs(«)Sv(r)dv, (27) 
Jo Jo 
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where c is t he speed of light, an d By (T) the Planck function at tem- 
perature r. Li & Draine (2001) show that when grains are bathed 



in the local ISRF of the Solar Neighbourhood, this continuum ap- 
proximation is valid for grain radii greater than 250 A. We adopt 
this same transition radius for the range of models considered here. 

The emissivities of these grains are calculated by assuming 
them to radiate as modified blackbodies with a single characteristic 
temperature: 

Ev = Cv,abs5v(7"). (28) 



7 EVOLUTION OF GAS MASS AND METALLICITY 

To follow the evolution of a star-forming system we use a simple 
closed box model. We denote the total stellar mass-loss rate (due 
to stellar winds and supernovae) as with metallicity Zoutflo„. 
The total mass of gas and dust in a cloud is nic(/), of which metals 
comprise a mass m^{t), such that Z(/) = m^{t) / m^{t) . Thus we may 
write: 



dmc 
"dT 



(29) 



We consider two models for the mixing of metal-enriched ma- 
terial in the cloud. In the simpler, we assume a time-invariant metal- 
licity, such that mz(0 = Zm(.ioud(0- We call this model the 'constant 
Z' mixing model. In the second, we assume perfect mixing, such 
that the metallicity of the material which goes into forming new 
stars is representative of that of the whole cloud, in which case: 



draz 
"dT 



-Z(f)1'(0+Zoutflow(f)2(f)- 



(30) 



We call this the 'perfect mixing' model. 

To calculate Z(f) and Zoutflow(0 used the same STAR- 
BURST99 models as described in Section^ The evolving metallic- 
ity of the gas forming into new stars could not be treated smoothly, 
as the Padova stellar-evolution tracks used by STARBURST99 are 
only available for five stellar metallicities: Z = 0.0004, 0.004, 
0.008, 0.02 and 0.05. Instead, we modelled each SSP within our 
calculation of the SFH using the track whose metallicity was loga- 
rithmically closest to that desired. 

As a simple starting point, we model two dwarf-galaxy- 
sized clouds. Firstly, we model a system composed initially of 
5 X 10^ Mq of gas and dust of metallicity 0.05, but with no stars. 
This cloud begins to undergo star formation at a rate of 0.3 yr^ ' 
at a redshiftn of z = 10, and maintains this rate of star formation 
until the current epoch. We show the evolution of the mass and 
metallicity of this system in Figure ^. 

Secondly, we form a model of an early-type galaxy. We as- 
sume a system of total mass 5 x 10^ Mq, composed initially of 
zero-metallicity gas. We assume it to undergo star formation at a 
constant rate in the redshift range 10-5, such that 90 per cent of 
this material is converted into stars over this period. The evolution 
of the mass and metallicity of this system is shown in Figure |5| 
We note that our perfect mixing model predicts such galaxies to 
have decreasing metallicities with time (Figure ^). After z ~ 5, 
mass loss from low-mass stars dominate the return of material to 
the ISM. This material has lower metallicity than supernova ejecta. 



See Section |l| for a definition of our adopted cosmology. 
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Figure 4. A model of the ISM of a galaxy of total mass 5 x lO' Mjj , initially 
of metallicity 0.05 solar, but containing no stars. It begins to form stars at a 
rate of 0.3 yr^' at z = 10, and maintains this rate of star formation until 
the current epoch, is the total mass of gas and dust in the galaxy's ISM. 
Panels (a) and (b) represent differing models of the mixing of material in 
the ISM, as described in the text. 



8 RESULTS : VARIATIONS IN DUST EMISSION WITH 
PHYSICAL CONDITIONS 

In this section, we present model results using the diffuse geometry 
described in Section 



3.3 



The simple optically-thin approximation 
allows us to investigate how the emissivity of dust depends upon the 
physical properties of the grain population, and the radiation field 
which heats the dust. Throughout this section, we take for our inter- 
stellar heating radiation field a STARBURST99 model of a continu- 
ous SFH of age 1 Gyr, and assume the ISM to be contained within a 
homogeneous spherical volume of radius rg. This leaves three free 
parameters: n^, x and vq. We choose, however, to re-parameterise 
the problem in terms of more physically interesting quantities - the 
column density of hydrogen nuclei along a diameter of the ISM, 
A'c, the total star-formation rate of the heating radiation field, "P, 
and - by means of the relations: 



'PL 



v,0 



(31) 
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Figure 5. A model of the ISM of an early-type galaxy of total mass 5 x 
10' M0, which converts 90 per cent of this material into stars in the redshift 
range 10-5. is the total mass of gas and dust in the galaxy's ISM and 
the mass of metals. Panels (a) and (b) represent differing models of the 
mixing of material in the ISM, as described in the text. 



and 



'iH =A'c/2ro, 



(32) 



where Ly o is the specific luminosity of a STARBURST99 model of 
unit SFR. The former relation may be derived by setting the second 
term on the right-hand side of Equation ( |l3[ ) to equal "PLy o- 

We truncate our adopted ISRF at an upper energy bound £max, 
which we set to the ionisation energy of hydrogen, 13.6eV; as- 
trophysically, we would expect inter-stellar neutral hydrogen to 
severely attenuate the radiation field at shorter wavelengths. To 
avoid altering the bolometric luminosity of the ISRF as a result 
of this truncation, we assume that the truncated radiation is repro- 
cessed to longer wavelengths with a spectrum matching that of the 
nebula continuum component of the STARBURST99 model. 

In Figure^ we illustrate the contributions made by grains com- 
posed of silicates, graphite, and PAH molecules to the total dust 
emission. The division between the contributions made by graphitic 
and PAH grains is somewhat blurred because of their being incor- 
porated into a single carbonaceous grain population with smoothly- 
varying optical properties. For the purposes of the decomposition 
shown in Figure W, we take 'graphitic' grains to be those with radii 



Total emission 

PAH emission 

Graphite emission 
- - Silicate emission 
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Figure 6. Model spectrum of the emission from dust in our diffuse geome- 
try, neglecting stellar emission. We assume an ISRF appropriate for ;f = 1, 
= lO^** Hm"', 'i' = 1 M, . , yr" ' , for a system of age 1 Gyr. We decompose 
the total emission into the contributions from grains composed of graphite, 
PAHs and silicates. 




\jixm 

Figure 7. Spectra modelled using our diffuse geometry with = 
102-*H m - ; we plot the sum of stellar and dust emission. We illustrate the 
effect of varying the intensity % of the ISRF with respect to that of the solar 
neighbourhood. 



> 50 A, and 'PAH molecules' to be those grains with radii < 50 A. 
We see that the emission shortwards of 20p.m is dominated by 
grains in the carbonaceous log-normal peaks, whilst the silicate 
grain population contributes around 70 per cent of the FIR thermal 
emission. The graphitic grain population is seen to play a relatively 
minor roleiii^sha^ing^^he^EDj^iiagr^^ 

& Draine ( p001| ) that the upper limit on the sizes of the carbona- 
ceous grains in their model was poorly constrained. 

In Figure ^ we show the effect of varying the intensity % of 
the ISRF with respect to that of the Solar Neighbourhood, whilst 
holding A'c and the star-formation rate "P of the galaxy constant. 
This means effectively altering the volume occupied by the dust 
and hence the energy density in the ISRF at fixed total luminos- 
ity. As % decreases, the peak of the FIR thermal emission is seen 
to move to longer wavelengths. The temperature of the grains con- 
tributing to this emission can be estimated from peak of the grey 
body emission, showing that the temperature of the large grains 
falls from about 24K to 13K as % varies from 10 to 0.2. 

In the MIR, the luminosities of the PAH features are essen- 
tially independent of %. This behaviour is explained by consider- 
ing the transient nature of the heating of the small grains. For the 
range of conditions considered, these small grains cool efficiently 
between photon absorption events. Thus, the MIR spectra of galax- 
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Figure 8. Spectra modelled using our diffuse geometry with Nc = 
10^'' Hm^^, including both stellar and dust emission. The effect of vary- 
ing the ionisation fraction / from its default of 0.8 is illustrated. 



ies are independent of whether a small number of grains are each 
undergoing frequent excitations, or a larger number of grains are 
each undergoing less frequent excitations. If we consider a single 
photon of stellar origin, the probability of its being absorbed de- 
pends only upon the integrated column density of dust it traverses 
in its passage through the ISM. Thus, the emission of transiently- 
heated grains depends only upon and the form of the spectral 
energy distribution and hence "P, independent of X- 

In Figure |8| we show the effect of changing the ionisation 
fraction / of the PAH grains, as defined in Section ^ from its de- 
fault of 80 per cent. For a discussion of observational evidence that 
PAH ionisation s tates vary between differing astrophysical environ- 
ments, see, e.g., Draine & Li (2001 ). We see that the MIR features 
in the 5 < A/ \im < 9 region are enhanced in ionised PAHs; those in 
the 10 < X/ \Lm < 20 region are enhanced when PAHs are neutral. 

The effect of PAH ionisation can be understood in terms of 
the adopted absorption cross sections shown in Figure ^ The en- 
hancement of the luminosity of the 5 < X/ [im < 9 PAH features 
when PAHs are ionised can be explained simply from our adop- 
tion of a C^^g («) which is larger at these wavelengths for ionised 
PAHs; the enhancement of the 3.3 [im feature when PAHs are neu- 
tral may similarly be explained. The enhancement of the features in 
the 10 < X/[im < 20 region when PAHs are neutral, however, can- 



not be so explained; C^^^ (a) is essentially independent of PAH 
ionisation state at these wavelengths. The explanation instead lies 
in the energetics of these grains. The UV-visible absorption cross 
sections of PAH molecules are essentially independent of their ion- 
isation states; neutral and ionised grains re-process UV-visible ra- 
diation into the MIR at essentially identical rates. If neutral PAHs 
are less luminous than ionised PAHs in the 5 < X/\im < 9 region, 
they must be more luminous at other wavelengths to have the same 
bolometric luminosities. 

We now examine the variation in the MIR feature strengths 
as a function of the physical properties of the population of PAH 
molecules. The size distribution of PAH grains may depend criti- 
cally upon environment; in H II regions, for example, strong fluxes 
of ionising photons may destroy many of the smallest grains. Such 
speculati on is s upported by the observation (see, e.g., Weingartner 
& Draine 200 1[ ) that a diverse set of grain size distributions are re- 
quired to fit Milky Way sight lines. 

In Figure l^al) we decrease the number be of interstellar car- 
bon nuclei per hydrogen nucleus locked up in PAH molecules from 
its default value of 6^ = 6 x 10^^ . We meanwhile hold the power- 
law component of the grain size distribution unchanged. We see 



that the result of removing some of the smallest grains from the 
model in this way is that the MIR luminosity decreases - eventually 
to leave only the Jeans tail of the stellar emission in its place. This 
result is as might be expected as carbonaceous grains with a < 15 A 
contribute in excess of 90 per cent of emission at A < 10 \im. 

In Figure ^bl), the effect of altering the relative proportions 
of the numbers of carbon nuclei composing PAH molecules in the 
two peaks is illustrated; as defined in Section ^ Z?c,i is the propor- 
tion in the 3.5 A peak and bc.2 = 1 — ^c,l that in the 30A peak. Our 
default model has bc.i = 0.75. We see that the effect this change on 
dust emission is to enhance emission at A < 17 ^.m when a greater 
fraction of the carbon nuclei are placed in the 3.5 A peak, while re- 
ducing emission at A > 17)i.m. A pivoting motion is seen around 
A ~ 17 )i.m, where the emissivities of grains in the two log-normal 
peaks, normalised per unit carbon atom, approximately equal one 
another. 

In Figure ^(cl), we consider the effect of varying - the 
grain radius, as defined in Equation ([T^), at which the optical prop- 
erties of the carbonaceous grain population shift from those of 
graphite to those of PAH molecules. We see that the features at 
A < 1 1 nm are little altered when is changed from 50 A to 10 A, 
but that further reducing to 4 A does suppress those PAH fea- 
tures with respect to the MIR continuum. As shown in Figure ^al) 
above, emission at A < 17 ^m is dominated by dust in the 3.5 A 
log-normal distribution, and so is only affected once is reduced 
to a comparable size. We conclude that our model pr edictions are 
relatively insensitive to , a result in agreement with Li & Draine 
( 2001 ), who found it to be relatively poorly constrained by obser- 
vation. 

Thus far w e have used as our default g rain size distribution 
that preferred by Weingartner & Draine ( 2001 ) for Milky Way sight 
lines with Ry = 3.1. These are typically sight lines which pass 
through infrared cirrus, away from dense molecular clouds. How- 
ever, those authors also fit grain size distributions to the dust along 
sight lines which pass through more dense parts of the Milky Way, 
characterised by larger values of Ry. For each value of Ry, they 
present a range of fits for different assumed values of fee- In Fig- 
ure [l^, we show the spectra produced by our model when we use 
the grain size distributions given by those authors for Ry = 4.0 and 
5.5, in each case using the value of be which they find to best fit 
the sight lines studied. For these enhanced values of Ry, we find 
that the infrared emission is reduced in line with what would be 
expected from the reduced number of dust grains per unit column 
density in these size distributions. 



9 RESULTS: CIRCUMNUCLEAR SHELLS 

In this section, we present spectra modelled using our circumnu- 
clear geometry. As in the previous section, our default heating ra- 
diation field is a stellar population with a continuous SFH, star- 
formation rate I* = 1 Mp, yr^ ' and age 1 Gyr. We again truncate 
this at an upper energy bound of £max = 13.6eV and scale up the 
nebula continuum component of the heating radiation field accord- 
ingly. 

We find that the predictions of our model depend only weakly 
upon the geometrical thickness, ri — ry, of the dust shell. For ex- 
ample, the emission from a dust shell of inner radius = 1 kpc 
and outer radius ri = 3 kpc is little different from that of an in- 
finitesimally thin shell at a radius of 2 kpc. Within the accuracy of 
our numerical integration, these two models are indistinguishable 
at A < 14 nm. At longer wavelengths, they do differ, the thicker 
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Figure 9. In each row, the left plot shows spectra modelled using our diffuse geometry with Nc = lO^'^Hm^^, X = ^ "I* = 1 Mq yr^'. These include both 
stellar and dust emission, and illustrate the effect of varying the the size distribution of carbonaceous grains upon the overall dust emission. The right plot 
shows the size distribution of carbonaceous grains used in the plot to its left. Panel (al) illustrates the effect of varying the strengths of the log-normal peaks 
in the size distribution. Panel (bl) illustrates the effect of varying the relative strengths of the two log-normal peaks. Panel (cl) shows the effect of varying 
the grain radius Og at which the carbonaceous grain population makes its smooth transition from having the optical properties of PAH molecules to those of 
macroscopic graphite particles. 




Figure 10. Spectra modelled using our diffuse geometry with = 
lifi'* Hm^^; we plot the sum of stellar and dust emis s ion. W e use the grain 



size distributions preferred by Weingartner & Draine ( 2001 ) for Milky Way 



sight lines with Ry = 3A, Ry = 4.0 and Ry = 5.5. The values of be in these 
models are, respectively, 6 x 10^^, 4 x 10"' and 3 x 10"'. 



shell having an enhanced luminosity at 15-70 (im, enhanced by 
23, 72 and 29 per cent at 20, 40 and 60|i.m respectively. This can 
be understood in terms of two effects: an enhanced probability of 
multiple-photon heating of small grains close to the inner edge of 
the shell, and radial variations in the temperatures of the large dust 
grains. For the remainder of this section, we adopt a geometry with 
ro = 1 kpc and ri = 3kpc. 

Figure [h] shows model spectra for dust shells with column 
densities A'c = 10^^ and 10^^ Hm"^; we assume solar metallicity 
for both the stellar populations and the dust shells. For each, we 
show the contribution made by attenuated stellar radiation - i.e. 



(1) 



in Equation I 



- given by: 



'^stellar 



: Lye 



(33) 



where Ly is the luminosity of the central heating source. 

When Nc = lO^'^Hm"^ stellar emission contributes 95 per 
cent of emission at 4^im. Convolution of this spectrum with the 
passbands of Spitzer's InfraRed Array Camera (IRAC) reveals that 
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Figure 11. Spectra of dusty star-forming regions modelled using our cir- 
cumnuclear geometry. For the heating radiation field, a stellar population 
with continuous SFH, star-formation rate 4* = 1 M0yr^' and age IGyr is 
used. Panels (a) and (b) show spectra of dust shells of column densities 
Nc = lO^'* and 10^* m^^ respectively. 



stellar emission contributes 92, 91, 34 and 1 1 per cent in the chan- 
nels at 3.5, 4.5 5.8 and 8.0 nm respectively. 

At Nc = 10^* Hm^^, these contributions reduce to 60 per cent 
of emission at 4 ^m, and 48, 45, 3.9 and 0.88 per cent in the IRAC 
channels at 3.5, 4.5, 5.8 and 8.0 [i-m respectively. An absorption 
feature in the attenuated stellar emission component is seen around 
10 Jim, indicating that the silicate 9.7ii.m feature has an optical 
depth of 0.4 in this model. 

In these model spectra, two PAH features appear well-isolated 
from their neighbours: that at 3.3 |i.m and the new 1.05 vim feature 
added by DL07. In the A'^ = 10^^ Hm^^ model, these have equiv- 
alent widths 0.3 and 0.04 \im. respectively. In the A'c = 10^^ Hm^^ 
model, the 3.3 \im. feature has equivalent width 0.07 \im; the 
1.05 \im feature is not apparent. 

In these plots we have, as previously , assum ed the grain size 
distribution which [Weingartner & Draine (2001) fit to Milky Way 
sight lines with Ry = 3.1, assuming foe = 6 x 10^^. From Fig- 
ure |lo| we can see that the effect of using a size distribution ap- 
propriate for a larger value of Ry would be a small reduction in 
dust emission at all wavelengths. 

In Figure |l^ we trace the evolution of our model spectra as a 
function of age for a range of SFHs. Figure |I^a) traces the evo- 
lution of the early-type galaxy modelled in Figure ^ Figure p^b) 
traces that of the system with ongoing star formation modelled in 
Figure ^. 

Figure |l2|(c) shows a system with the same SFH as that mod- 
elled in Figure [l2(b), but assuming a constant dust shell column 
density A'c = lO^Hm^^. The time-invariance of the emission at 
X < 300 nm can be understood in terms of an equilibrium being 
reached between star formation and stellar death among the mas- 
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Figure 12. The time-evolution of the spectra of dust-enshrouded star- 
forming regions modelled using our circumnuclear geometry. In all cases, 
ro = 1 kpc and r\ = 3 kpc. Panel (a) shows the evolution of the SED of the 
early-type galaxy modelled in Figure ^ the star-formation history and dust 
column density assumed are as shown in that figure for the perfect mixing 
model. Panel (b) shows the evolution of the system with ongoing star for- 
mation modelled in Figure ^. Panel (c) shows a system with the same SFH 
as Panel (b), but assuming a constant column density = 10^' Hm"^. 



sive stars which produce it. This contrasts with the behaviour seen 
in Figure |l^b), where the system becomes more metal-rich and 
dusty as it evolves, increasing extinction at ultraviolet wavelengths. 



10 THE DEPENDENCE OF MIR LUMINOSITY UPON 
SFR 

To probe the dependence of the MIR luminosities of our models 
upon SFR, we constructed a set of circumnuclear models, each hav- 
ing a continuous SFH of age 1 Gyr but different star formation rates, 
ranging from 0.02 to lOOMi^yr^'. We adopt otherwise constant 
model parameters: ro = Ikpc, r\ = 3 kpc and A'c = 10^^ Hm^^. We 
constructed two sets of models, using the LDOl and the DL07 PAH 
cross sections, to compare their predictions. 

We compare these to the sources in the Spitzer extragalac- 
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tic First Look Survey (FLS) field. We correlated this source cat- 
alogue with the Fourth Data Release of the Sloan Digital Sky Sur- 
vey (SDSS;|York et alj|200Cl|; [Adelman-McCarthy et al]|200^) us 



ing a 2" matching radius to obtain redshifts for these sources, and 
used emission line measurements from the MPA/JHU catalogue^ 
To remove AGN and spurious/marginal detections from the cata- 
logue, and to correc t for a perture effects and extinction, we fol- 
(2004): AGN were identified u si ng th e spec- 



lowed Nikolic et al. 



tral classification scheme of Veilleux & Osterbrock (1987); aper- 



ture and extinction correctio ns were applied to th e emission line 
data using the prescription of Hopkins et al, ( 2003 ). We estimated 
the SFRs of the sources by comparing their Ha luminosities with 
that predicted by STARBURST99 for a 1 Gyr-old continuous SFHs. 
A'-corrections for the Spitzer 8 ii.m and 24 \im fluxes were derived 
using our models described above, selecting for each object that 
model which most closely matched the Ha-derived SFR. 

We convolved our model spectra with the 8-|J.m and 24-\im 
passbands of Spitzer to produce matching data from our model and 
the observed source population. The results are shown in Figure |l3[ 
also shown are power-law fits to the observed source population, 
the fit coefficients for which are listed in Table ^ (c.f. a similar fit 
performed by Wu et al 2005 ). In the upper panels, we use the LDOl 



PAH cross sections (see Section in our models; in the lower pan- 
els we use the DL07 cross sections. The non-linearity of our model 
Ly {24 ]im) — relations at > lOMQyr^' results from the in- 
creasing emission of grains excited by multiple photons at the high- 
est star formation rates. 

These results are insensitive to the chosen geometry, as the 
emission from our models at both 8 and 24 ^.m is dominated by 
grains in the transiently heated regime. The luminosities are, how- 
ever, dependent upon the column density of dust; this is likely to 
be one cause of the scatter of the FLS sources from our model re- 
lations. Increasing Nc by an order of magnitude, to 10^^ Hm^^, in- 
creases our predictions of Lv(8 \xm) by a factor of around 1.75, and 
Lv(24 p.m) by a factor of around 2. Increasing further causes our 
predictions for Lv(8vim) to decrease as self-absorption becomes 
significant. Further scatter in these plots is likely to result from the 
uncertainly in the conversion of Ha luminosities into star forma- 
tion rates. 

It is interesting to note that a better fit to the 24-fi.m lumi- 
nosities of the FLS sources is obtained using the LDOI PAH cross 
sections compared to the DL07 cross sections (see Figures [l3| ? and 
[l3|d). This difference can be explained by the enhanced PAH cross 
sections at 24jim in the LDOI model, as seen in Fig ures ^(a) and 
^b). The fit to the 24 [xm luminosities in Figure |l3|(d) can be im- 
proved by increasing Nc, but this leads to an over- prediction of the 
8-^im luminosities in Figure pj|(c). 
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9.45 ±0.05 


0.87 ±0.08 
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9.06 ±0.04 


0.97 ±0.07 
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Table 2. Coefficients of the power-law fits shown in Figure where 
log 10 (vLv(A)/Lr^i) = a±fe X logiQ (l'/M,.,yr"'). The number of sources 
used in each fit is N. Although the A^-corrections derived for the FLS galax- 
ies differ depending upon whether the LDOI or DL07 PAH cross sections 
are used, these differences are modest; the changes to the power-law coeffi- 
cients above are much smaller than the quoted en'ors. 
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Figure 14. Plot of the positions of galaxies modelled by our circumnuclear 
geometry in colour-colour space. Each trace represents the evolution of a 
model galaxy from an age of 1.6Gyr (top-right end) to 13.8Gyr (lower- 
left end). Late-type galaxies, with continuous SFHs, are shown with solid 
traces; early-type galaxies, with instantaneous SFHs, are shown with dashed 
traces. Three traces are shown for each galaxy type, representing galaxies 
enshi'ouded with three different column densities of dust. All colours ai'e 
shown in magnitudes, and are calculated in the galaxy rest frame. 



11 PREDICTED GALAXY COLOURS 

Using the models presented in Section ^ we can make predic- 
tions of the relative merits of colour diagnostics in the discrimina- 
tion of systems with ongoing star formation from those with only 
old stellar populations. In Figure |l^ we show the colours of late- 
and early-type models, with continuous and instantaneous SFHs re- 
spectively, enshrouded by three different column densities of dust. 
In the lower panel we consider colours formed between the /-band 
of the UK InfraRed Telescope (UKIRT), as used by the UKIRT 



Available from: ittp://www. mpa-garching.mpg.de/SDSS/ 



Infrared Deep Sky Survey (UKIDSS), in conjunction with Spitzer 
MIR photometry. In the upper panel, we consider colours which 
use only IRAC photometry. In all cases, the late-type models ap- 
pear redder than systems with only old stellar populations, a result 
of the enhanced UV- pumping of PAH grains. Models which are en- 
shrouded by higher column densities of dust appear systematically 
redder, as a result of their enhanced dust emission. 

All of these colour predictions are robust to changes in star- 
formation rate *P. In our model, both stellar and PAH emission are 
linearly proportional to 1* except at the highest star formation rates 
(see Figure O); only the far-infrared greybody varies in colour with 
"P. These predictions are also robust to changes in the assumed 
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Figure 13. Solid lines: The Spteer-convolved 8-mn (Left panels) and 24-|am (Right panels) luminosities of circumnuclear models with continuous SFHs, 
age IGyr, ro = Ikpc, ri = 3kpc and Wc = lO^^Hm^^. Dashed lines: Power-law fits to the Spitzer extragalactic FLS source population (points). In the upper 
panels we use the LDOl PAH cross sections; in the lower panels we use the DL07 cross sections. 



masses of the galaxies if is assumed to scale linearly with mass 
whilst remains constant. In Figure |l5| we investigate the red- 
shift and column density dependence of the colours shown in Fig- 
ure by tracing the redshift evolution of systems of constant N^^. 
The diagnostic power of colours formed between IRAC channels 
is limited beyond z > 0.5; at higher redshifts these bands are in- 
creasingly dominated by stellar emission, especially at lower dust 
column densities. The power of 7 — 8 ii.m and J — 24 [im colours 
extends to higher redshifts, though a degeneracy exists between 
highly obscured old stellar populations and recent star formation. 
In all cases, the two SFHs considered in Figure |l^ converge at high 
redshift; both contain relatively young stellar populations shortly 
after their formation at z = 10. 



12 SUMMARY 

We have developed a semi-empirical model of the emission of dust 
in star-forming galaxies, capable of producing model spectra which 
extend from the Lyman limit through three orders of magnitude in 
wavelength to the far-infrared. A particular strength is its state-of- 
the-art model of the mid-infrared PAH features: we treat the tran- 
sient heating of sma ll dust grains using the sophisticated method of 
Draine & Li|(|200l|) and use the grain interaction cross sections of 



Draine & Li^ 



In Section H, we showed how the predictions of our model de- 
pend upon the microscopic properties of the dust grain population 
- its grain size spectrum and ionisation state, for example. In Sec- 
tions ^ and |l^ we went on to demonstrate the spectra produced by 
our model depend upon the large scale properties of galaxies: their 
star formation histories and dust column densities. In Figure |l3[ we 
found that the Spitzer extragalactic First Look Survey (FLS) source 
population is better matched by models which use the PAH cross 



sections of Li & Draine teooi| ) than those which use the cross sec- 
tions of praine &Li| ( [Z007| ). 

In Section |ll|, we provide estimates of the near-to-mid- 
infrared colours of early- and late-type galaxies, demonstrating the 
power of our model as a tool for interpreting colour-colour dia- 
grams. 
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APPENDIX A: CHOICE OF ENERGY BINS 



The energy bins used in Section 5.1 are chosen so that P{E) is 



well-sampled in regions where it is changing rapidly. For the small- 
est grains, this means that the ground state, and those close to it, 
must be well-sampled. To ensure that this is the case, we set the A'l 
lowest-energy bins in our scheme to contain the Ni lowest-energy 
vibrational states of dust grains, one in each bin, thus treating them 
exactly. At the opposite extreme, we require the occupation prob- 
ability of our highest-energy bin to be very small, such that the 
high-energy tail of P{E) is not truncated. 

To meet these requirements, we place the next A'2 bins at linear 
intervals of Un^ —Uni-i above U ai, . We then place the remaining 
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Figure 15. The colours of galaxies modelled by our circumnucleai' geometry with continuous (left panels) and instantaneous (right panels) SFHs; in each case, 
star formation starts at z = 10. Each trace represents the redshift evolution of a system of constant N^. All colours are shown in magnitudes, and are calculated 
in the frame of the observer. 



^bins ^ N2 bins at logarithmic intervals above Uf^^^f^-, with 
multiplicative spacing a given by: 



In this paper, we use A'^i = 20, A'2 = 50 and A'^i, 



600. 
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APPENDIX B: 
MATRIX T 



CALCULATING THE TRANSITION 



Several approximations for calculating the emission from stochaS' 



t ically heated dust grains are discussed in detail by Draine & Li 



(2001). They all rely on dividing the possible enthalpy content of 
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the grain into of the order of 500 bins and calculating the probabil- 
ity of a grain transitioning from an enthalpy that belongs to one bin 
to an enthalpy belonging to a different bin. 



2001 



The upward transition rates Tji are given by (Draine & Li 
equation 15): 



^J' = 77 TT / G;,-(£)Cabs (E)uE {E) dE for ; < N^ns, 

U j U i JWi 



M V 



E-Wi 
Wc-Wi 



+ / C.^bs{E)uE{E)dE 



Qbs(£)"£(£)d£ 



for j = A^bins 



(Bl) 



(B2) 



where Qbs is the grain absorption cross section, ue is the energy 
density of the radiation field and Gj, is a correction factor for the 
finite width of bins and is given by: 



AUjAUi 



Gji = { 



min {AUj,AUi) 
AUjAUi 
W4-E 



AUjAUi 







Wi <E < W2, 
W2<E<W2, 

W2<E< W4, 
otherwise. 



(B3) 



The integration limit quantities Wi, W2, W3 and W4 are defined in 
Figure gl|. is equal to Uj^^^^^^ - C/™". 

Throughout this paper, we use the c ontinuous cooling a pprox- 
imation for the downward transitions (Draine & Li 2001, equa- 
tion 41): 

T _ 1 8;r 



Ui-Uj h^c^ 

E^C^^^E) 



X / Eu 



exp{E/k6u) - 1 



dE 



fori> 1;7 = !-1(B4) 



T;V = 



for all other {i,j} 



(B5) 



where 0„ is t he ch aracteristic temperature of bin u, as defined by 



Draine & Li (2001) 



The computational cost of evaluating the transition matrix Tji 
is dominated by evaluating the upward transitions, in particular the 
integrand in Equation (Bll), separately for each {/, j} pair. It may 
be reduced by noting that, although G^,- is a function of ; and j, 
it always remains a linear function of E. Hence the integrand in 
Equation (|bi|) is, for any {i,j} pair, a linear combination of: 



/] {a,b) ■■ 



rb 

/ dEC.^bs{E)uE{E) 

■la 



(B6) 



and 



h{a,b) 



f''dE-E-C,u{E)uE(E). (B7) 

Ja 



where Ii and I2 are independent of / and j. We can then write 
Tji as a linear combination of h {Wi , W2), h (^2,^3), h (^^5,^4), 
h {W\ , ) I ^2 ( W2 , 1^3 ) and /2 (1^3 , VV4 ) , with coefficients which are 
independent of E. We evaluate /[ and I2 by non-adaptive integration 
using a pre-computed grid of the two integrands appearing in 1\ and 
h- By using this grid, we reduce the number of necessary evalua- 
tions of Cabs(£') and ue{E) from ff{N^) to ^(10^). Non-adaptive 
integration is appropriate since ue, the output of stellar population 
models, is computed on a fixed grid. 



W2 
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^max 



jymin 
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Figure Bl. Illustration of the energy boundaries used in function Gji{E), 
where and (' are the indices of the upper and lower bins respectively and 
J/™'" and are the boundaries of the rth bin. Note that W2 is the smaller 
and W3 the larger of the pair of values: - J7™) and ({/)"'" - Uf"). 
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